Seebeck effect in the graphene-superconductor junction 



Marcin M. Wysokiriski^B and Jozef Spalek 1 - 2 >Q 

1 Marian Smoluchowski Institute of Physics, Jagiellonian University, Reymonta 4, PL-30-059 Krakow, Poland 

^Faculty of Physics and Applied Computer Science, 
AGH University of Science and Technology, Reymonta 19, PL-30-059 Krakow, Poland 

(Dated: February 5, 2013) 

Thermopower of graphene-superconductor (GS) junction is analyzed within the extended Blonder- 
Tinkham-Klapwijk formalism. Within this approach we have also calculated the temperature de- 
pendence of the zero-bias conductance for GS junction. Both quantities reflect quasi-relativistic 
nature of massless Dirac fermions in graphene. Both, the linear and the non-linear regimes are 
considered. For the sake of completeness, in Appendix [B] we discuss limitations and consequences 
resulting from different approximations in the transverse-mode counting. 



I. INTRODUCTION 

Graphene is one of the most remarkable new materi- 
als. Not only has its discovery^ violated in same sense 
Landau's theory of the thermodynamical instability of 
a two-dimensional structure^, but also due to the pecu- 
liar band structure it has provided us with an invaluable 
opportunity to test relativistic quantum electrodynamics 
in the desktop laboratory 3 -. For that reason, much ef- 
fort has been put into understanding of the phenomena 
associated with this material. 

Since the graphene-based devices are usually consid- 
ered as a mcsoscopic systems, the Landauer approach is 
widely utilized to study the ballistic transport in them^r— . 
Even though this approach does not account in the sim- 
plest form for all the features of the material, it pro- 
vides a good overall description of the electric trans- 
port. This approach was extended by Blonder, Tin- 
kham, and Klapwijki (BTK) to the case of standard 
normal metal-superconductor junction (NS), and in this 
manner yielded a very good description of experimental 
data^. Their method has been widely used for different 
specific situations^— and finally adapted for graphene- 
superconductor hybrid systems^—. One of the most 
peculiar properties predicted in such systems is the spec- 
ular Andreev reflection and deviations of the conduc- 
tance spectra— from those predicted by BTK for normal 
metals^. 

Landauer formalism has also been successfully adapted 
for thermoelectrical transport in mesoscopic devices^—. 
In the case of standard NS junctions BTK formula, 
also turned out to be useful technique for predicting ef- 
fects concerning thermal properties of electric and heat 
current o 11 ' 13 ' 17 ! 26 . This method has also been used for 
the graphene-based superconducting hybrid structures 
for obtaining the thermal conductanc e 20 ' 27 " — . However, 
the thermopower has not been studied so far. This topic 
is addressed in this article. 

In this work we provide systematic study of the effect 
of the temperature on the charge current in the super- 
conducting graphene junction (GS) using a generalized 
BTK formalism for the specific case of graphene. We 
present results concerning the temperature dependence 



of the zero-bias conductance and the Seebeck coefficient 
in the linear regime. For the sake of completeness, we 
also discuss the non-linear thermopower. 

The paper is organized as follows. In Sec. [TT| (and in 
Appendix [X]) , we present briefly a generalized BTK ap- 
proach for the charge current through the GS junction. 
The linear transport coefficients are discussed, in partic- 
ular the zero-bias conductance and the thermopower. We 
also briefly comment on the effect of non-linear correc- 
tions on the Seebeck coefficient. Finally, we conclude in 
the Sec. HVl As a supplement, in Appendix [B] we add a 
brief discussion on the transverse-mode counting. 



II. MODEL 

We consider a ballistic limit for graphene based junc- 
tion composed of the normal region and induced by 
means of the proximity effect superconducting region (cf. 
FigO}. For the description of the unconventional quasi- 
particle states we utilize Dirac - Bogoliubov - de Gcnncs 
equations for the two-dimensional (2D) sheet of graphene 
in the formic 



Hj -E F 1 
A+ 



A 

E F 1 - Hj 



(1) 



where the index j can be either + or - what refers to the 
two incquivalcnt valleys K and K' in the Brillouin zone. 
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FIG. 1. (Color online) Proposed schematic, experimental 
setup considered in our modeling. 



2 



The single particle Hamiltonian is given by 

H± = -HVF(CT x d X ± (JyOy) + U, 



(2) 



where vp is the energy independent Fermi velocity for 
the graphene, and {(Ji\ denote respective Pauli matrices. 
Because of the valley degeneracy one can effectively do 
the calculation for the one valley only. We assume that 
in the geometry, where the interface is determined by 
the y-axis, the pair potential with the s-wave symmetry 
changes step-like in the x-axis direction, 



A(r,T) 





A{T)e l 



x < 0, 
x > 0, 



(3) 



where the temperature dependence of the gap function 
can be deduced from the usual BCS theory^ and is taken 
in the following form 



A(T) 



= tanh 



(4) 



The BCS theory of the superconductivity based on the 
requirement that the coherence length is large when com- 
pared to the Fermi wavelength. Under that condition we 
can assume that the additional potential U has the form, 



U(r) 





-U 



x < 0, 
x > 0, 



(5) 



with large Uq (Ep + Uq >> Ao) and for the simplicity 
we define E' F = Ef + Uq. In numerical calculations we 
set E' F = 1000A . 

In the spirit of the BTK scheme (by matching the 
wave functions at the boundary (x = 0)), we obtain 
expressions for the amplitudes of the Andreev hole re- 
flection (AR) (a(e, 8)) and the normal reflection (b(e, 9)) 
of the incident electron - see Appendix [A] for details of 
the method. Note that there is no intrinsic barrier at 
the GS junction, and thus the Fermi vector mismatch 
is the source of the normal reflection. The transmission 
probability, averaged over the angles, takes the following 
formic, 



TOO 



r/2 2 



l&M)| 2 + ^&a(e,0)| 



(6) 

The BTK formalism combined with the specific trans- 
mission probability derived for graphene, defines charge 
current through the GS interface as^i&, 

h = T j_Je N(e)T(e) (f G (e - eV) - f s (e)) , (7) 

where f G , f s are the Fermi distribution functions for 
the normal (G) and the superconducting (S) region of 
GS junction respectively, and 



N(e) 



\E F + e\W 
itHvf 



(8) 



is the energy dependent number of transverse modes in 
the graphene sheet of width W— . However, formula (O 
is not always accurate. The additional assumption is 
needed, that for each mode carrying incident electron, 
having energy Ep + e and being Andreev reflected there 
is always enough modes at the level Ef — £ for this process 
to happen. For a more detailed discussion of this topic 
see Appendix [Bl 

The quantity describing thermoelectric properties of 
the system is thermopower, or Seeback coefficient (S) 
measuring the voltage driving to zero the current flowing 
in response to the temperature difference, namely 



S 



V 
ST 



(9) 



III. CHARGE TRANSPORT 
A. Linear regime 

Expansion of the Fermi functions in the normal and 
superconducting regions to the first (linear) order in both 
the bias and the temperature difference, with the average 
temperature T, i.e., 



f G = h-8T/2(e - eV) ~ / r (e) - eV^- + ^e^, 



f S = It+st/2^) - Me) - 



df_ ST df_ 
de + 2T 6 de 
Se df_ 



enables to decouple Eq.© in the form, 
L = GV + iJST. 



(10) 
(11) 



The resulting from it the linear transport coefficients can 
be thus rewritten in the following closed forms, 



G 



4e2 
h 



de %N(e)T(e), 
de 



(12) 



The temperature gradient and the bias are set as positive 
with respect to x coordinate. For the temperature in the 
system approaching zero, the expression for the electric 
conductance (G) reduces to the well-known BTK zero- 
bias conductance formula^ 



Gt^o = ^-JV(0)T(0). 



(13) 



In Figl2] we have plotted the zero-bias differential con- 
ductance (calculated from Eq. (fT2)) ') as a function of tem- 
perature, normalized by the ballistic conductance go with 
having TV transverse modes in a sheet of graphene of 
width W and given by 



4e 2 f°° df 



(14) 
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FIG. 2. (Color online) Normalized zero bias conductance as 
a function of temperature for various Fermi vector mismatch. 



For relatively high position of the Fermi level in graphcne 
(roughly Ep > 5Ao), the influence of the increasing tem- 
perature (up to the approximately T/T c = 0.5) on the 
conductance is almost negligible (cf. Fig. [2}. The origin 
of this behavior is strictly connected with the relatively 
slow variation of the transmission probability T(e) for 
the subgap energies in this temperature regime, and thus 
does not differ qualitatively from the standard NS case^. 
In the low doping regime (Ep < 5Ao) the linear differen- 
tial conductance as a function of the temperature vastly 
differs from the standard NS case and reflects the specific 
electronic nature of graphene. 

The non-trivial behavior of the Dirac fcrmions in 
graphene has also its influence on the linear ther- 
mopower. From Eq. ([9]) and formula for this quan- 
tity reads 

q= lJ_ 1 Co.de e\e + Ep\%T(e) k B 

G k B T J^dele + Ep^Tie) e" 1 j 

Results obtained by the numerical integration are pre- 
sented in the Fig. [3l The low temperature regime dif- 
fers from the one obtained for the NS junction^. Con- 
trary to the NS case, for the graphene-based structure 
the thermopower does not vanish for non-zero temper- 
atures. This is due to the relativistic nature of charge 
carriers in graphene, where AR does not vanish even for 
the high effective barrier (Fermi velocity mismatch in our 
case) for the subgap energies (for the normal incidence 
AR happens always with certainty^). For the low dop- 
ing regime (Ep < Ao) we observe a clear maximum of 
the thermopower in the temperature range where the su- 
perconductivity exists. In the linear regime, the ther- 
mopower for the selected range of parameters near Tq 
is large and reaches values up to lks/e. This suggests 
that there is a potential for application of this setup for 
cooling of various nanostructures. 



FIG. 3. (Color online) Linear thermopower as a function of 
temperature in the system for various Fermi vector mismatch. 



B. Effect of the non-linearity 

We have also studied numerically the effect of the non- 
linearity in our system. The thermopower in the non- 
linear regime can be calculated as a ratio between the bias 
voltage and the temperature gradient, when no charge 
current is flowing (cf. Eq.©). The results in the non- 
linear regime are presented in Fig. |4l We have found 
that the non-linearity influences, the thermopower in not 
systematical manner with changing Fermi level position 
and the average temperature of the system. However 
the change is not as dramatic as in the NS case^ and is 
almost unnoticcablc in the doped regime Ep > 5A . 



IV. CONCLUSIONS 

In this work we have analyzed the thermoelectric 
charge transport in the graphene junction consisting of 
the normal and the superconducting parts. In the linear 
regime we have calculated the temperature dependence 
of the zero-bias conductance and the thermopower. We 
have found deviations of these quantities from the stan- 
dard normal metal-superconductor junction case that are 
caused by the relativistic nature of electrons in graphene. 

We have also studied the effect of non-linearity on the 
thermopower and we have found that for a high Fermi 
level positions (Ep > Aq) it stays almost unaffected and 
in the low Fermi level regime is noticeably enhanced with 
the increase of the temperature gradient. 

As a supplemental discussion in the Appendix [B] wc 
discuss some consequences of including into the standard 
BTK model the energy dependent number of transverse 
modes. We point out that in graphene with the Fermi 
level near the Dirac point, where the energy dependence 
of the number of transverse modes is not negligible, one 
needs to be cautious in using the BTK formula as ap- 
plied to the discussion of the effects associated with an 
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FIG. 4. (Color online) Nonlinear thermopower as a function 
of the temperature gradient set over the junction for vari- 
ous Fermi energies in graphene. The net temperature in the 
system is marked above each plot. 



interplay of the heat and the electric currents, since one 
of the Onsager reciprocity relation can be violated. 
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Appendix A: BTK for graphene 

The wave function in the normal part of graphene (NG) 
^jv and in the superconducting region (SG) ipg, look re- 
spectively as follows 



bip% +aijj h 



N 



(Al) 



where the superscripts e, h refers to electron and hole in 
NG and electronlike and holelike excitation in SG, and 
the superscripts + and - to right and left moving particle 
respectively. 

Spinors resulting from E q.CD ) are expressed in the sim- 
ilar manner as in the Ref. 1 1 QL i.e., in the form 



e± 
N 



[l,±e ±ie ,0,0] T e ±ifeexcos9 , 
[O,Q,l,e- l0A ] T e- lkhxcoseA , 



[u,ue s ,ve 
-ve s ,ue 



-i<t> V gi(S|-*)jT e iq e a:cos6l| 

10 , -ue~' l{e s + ^] T e~ iqhx cos0 s , 

(A2) 

where for the sake of clarity we do not include phase 
factor e yV since it corresponds to conservation of mo- 
mentum in y direction. The corresponding wave vectors 
are defined as follows, 



k e(h) = 



A 2 



Hv f Hvf 
and the coherence factors are given by 



(A3) 




Ve 2 - A 2 




The conservation of momentum at the interface and 
along y direction enables us to obtain mutual relations 
for the specific angles, namely 



k e sin (9 = k h sin* 



q e sin 6% = q h sin6>§. 



(A5) 



The system must also satisfy the continuity condition 
at the interface, Vvl(0) = Vv.r(0). The Hamiltonian is 
linear therefore is no need in matching derivatives. The 
resulting wave function amplitudes take the form, 
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Appendix B: On transverse-mode counting and 
fulfilling the Onsager relations 

In this Appendix, which has rather a Supplementary 
character, we raise the question, whether standard ways 
of including the number of transverse modes in the BTK 
formalism is applicable in the situation, when one consid- 
ers not only electric current, but also the heat transport 
as well. In particular, we address the question of the va- 
lidity of the Onsager relations which in our system should 
be obeyed^ 2 -. 



B.l. Two different attempts 

One encounters some difficulties when tries correctly 
considering number of transverse modes in BTK model 
for graphene. In the BTK schemed in order to obtain 
transmission probability one has to take the difference 
between distributions of the populations of particles mov- 
ing towards and away the interface, 

HO-/-.-/.- 

(1 -£)/-- Af+-(1- A- B)f, 

where for convenience and clarity we omit angle depen- 
dence and used generalized notation - A stands for prob- 
ability of the Andrcev reflection and B for the normal re- 
flection, and we define / ± = /(e±eV), whereas / = /(e). 
The above formula (|B1|) usually is used in the more com- 
pact form as 



T(e) = (l-B + A)(f--f). 



(B2) 



However, this is not a true transmission amplitude across 
the interface. Nevertheless, when one takes the integral 
over whole range of energies, it provides the same result 
as the true transmission 7~~'(e), 



/OO />00 
deT(e)= / deT'(e). 
-OO J — OO 



(B3) 



The above relation becomes not valid, when one in- 
cludes directly the number of transverse modes (which 
would not be an even function of e) under the integral. 
It is clearly seen when one calculates the resulting differ- 
ential conductance in the zero temperature for the trans- 
mission T(e), 

4e 2 d f°° 
G=^^\ deT(e)N(e) 

-oo (B4) 
= + A (eV) - B (eV))N(eV), 



h deV 



and analogically for the transmission 7~~'(e), 
4e 2 d 



G'=- 



h deV 



de T'(e)N(e) = 



4e 



(B5) 



= — [(1 - B(eV))N(eV) + A(eV)N(-eV)] . 
If N(e) is not an even function of e, G ^ G'. 



B.2. Standard approximation 

Usual way to include the number of transverse modes 
and utilized in majority of papers studying differential 
conductance of graphen o 18 ' 33 ' 34 is to take into account 
the transmission T(e) (Eq. IB2[) what leads to the formula 
(|B4[) . Physically it means that there is always enough 
modes for propagation of created in the AR hole beneath 
the Fermi level (for convenience, we call it Assumption 
1). The need of this assumption is clearly visible in the 
perfect AR regime. For ideal transparent interface and 
for subgap energies, AR happens with certainty. It means 
that in the normal region every incident electron should 
be converted into the reflected hole. However, due to 
mismatch in the number of available modes above and 
beneath the Fermi level, a part of incident electrons can- 
not Andrecv reflect. Since the transmission to the super- 
conductor is also forbidden, the excess electrons have to 
scatter via the normal-reflection process. 

In systems where the Fermi energy is quite high it does 
not matter for studying the electric current since the 
normal reflection due to the mismatch in modes num- 
ber above and beneath Fermi level would be negligible 
small. However, in the case of graphene in the low dop- 
ing regime, near the perfect AR, it can make a difference. 

The Assumption 1 is also used in many papers to study 
the thermal conductanc e 20 ' 27 " — . Therefore, a natural 
question arises, whether the Onsager reciprocity relations 
are satisfied (they should be obeyed in our system 3 - 2 -). 

For the ballistic transport and with the use of Assump- 
tion I, the heat current through the barrier is defined as 
follows, 



lQ= hJ OQ (/ G ( £ - eV ) / 5 ( £ )) ■ 

(B6) 

Note that here the transmission probability is different 
than in the electric current case, i.e., 



T"(e) = (1-B-A) (/"-/). 



(B7) 



Physically, it means that there is no heat current into 
the superconducting region for the subgap energies. The 
same form of the transmission (Eq. (|B7[) ) has been used 
for instance in Refs. @, [27rtI3 • 

In general we can write that the electric and the heat 
currents are driven by the bias voltage or the temperature 
gradient, so that 



I e = GV + iJST, 
la = IqV + KST. 



(B8) 



In that situation, the selected linear coefficients read, 



4e 



h 



df 



I Q=—T * -i-eN(e)T"{e) 



de 



(B9) 
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One can see clearly that within the Assumption 1 
Iq 7^ —Tf£, and thus the corresponding Onsager re- 
lation is broken. This is due to the difference between 
T" and T ■ Therefore, one should be careful in using 
this formalism for calculating for instance the thermal 
conductanc e) 20 ' 28 ' 29 , which from definition should results 
from an interplay between electric and heat currents, 
where the Onsager relations should be satisfied. 

B.3. Restoration of Onsager relations 

Contrary to the Assumption 1 one can take into ac- 
count full form of transmission (Eq. (|Bip ) and directly 
add the number of transverse modes under the integral 
to formulas for the electric and the heat currents. This 
attempt was used before in Ref. [Tl] and [2l|. Within this 
assumption one restores the Onsager reciprocity relation. 
Despite an apparent success, one pays a price of 'decou- 
pling' the coherent process of AR. More precisely, AR is a 
process in which incident electron in some sense borrows 
a partner electron from beneath the Fermi surface, and 
as a Cooper pair is absorbed into a superconductor at 
the Fermi level, and simultaneously the hole is reflected 
back into graphene. Therefore, it should be clear that 
there is the same number of incident electrons scattered 
in AR as holes reflected. In this attempt of including 



number of transverse mode, the reflection back is not 
guarantied. It is clearly visible if one considers the heat 
current through the interface in the subgap regime. Any 
energy transported towards the interface is taken back 
by reflected hole or electron, and therefore, the heat cur- 
rent should be zero. However, if we calculate within this 
attempt (with T") heat current, we obtain a nonzero en- 
ergy transported through the interface, what is clearly 
wrong. Note that still this way of including number of 
transverse modes is reasonable for the electric current for 
the high Fermi level position as in Refjlll. 

B.4. Remarks 

The Assumption 1 represents a reasonable approxima- 
tion for the electric-current calculation. Therefore, we 
use it for calculating the thermopower and the zero-bias 
conductance of the GS junction in this article. Further- 
more, the resulting formula for the charge flux is rigorous 
away from the perfect AR limit (where AR happen with 
certainty). 

However, since within this assumption the Onsager 
reciprocity relations in the linear regime are not fulfilled, 
for investigating the effects of intcrplaying the electric 
and the heat fluxes one should be careful to construct a 
model which would also obey Onsager relations. 
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